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Abstract 

We study the prospects of pinning down the effects of non-standard antineutrino interactions 
in the source and in the detector at the Daya Bay neutrino facility. It is well known that if the 
non-standard interactions in the detection process are of the same type as those in the production, 
their net effect can be subsumed into a mere shift in the measured value of the leptonic mixing 
angle ^13. Relaxing this assumption, the ratio of the antineutrino spectra measured by the Daya 
Bay far and near detectors is distorted in a characteristic way, and good fits based on the standard 
oscillation hypothesis are no longer viable. We show that, under certain conditions, three years of 
Daya Bay running can be sufficient to provide a clear hint of non-standard neutrino physics. 
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I. INTRODUCTION 



The past and ongoing neutrino oscillation experiments provide a firm evidence that the 
neutrino flavour is changing throughout the neutrino propagation. Except for some re- 
cent signals reported by MINOS jl| and MiniBooNE ^ (and also previously by LSND jsl) 
the vast majority of the data is consistent with the hypothesis of neutrino flavour oscil- 
lations driven by a pair of mass-squared differences: Amg^ = 2.45 ± 0.09 x 10"'^ eV^ 
(or Amg^ = — (2.34^QQg) x 10~^ eV^ if an inverse neutrino mass hierarchy is realized), 
often called the atmospheric mass-squared difference, and the solar mass-squared differ- 
ence Am2i = 7.59lo'^8 x 10~^ eV^, together with a pair of the corresponding mixing 
angles: sin^ 623 = 0.51 ± 0.06 (or sin^ 623 = 0.52 ± 0.06 the inverse hierarchy case) and 
sin^6'i2 = 0.312^q;qJ5 {4]. This, however, requires at least two of the oscillating neutrinos to 
be massive. 

By construction, neutrinos are massless in the Standard Model (SM). Thus, a lot of 
effort has been spent on devising SM extensions that could not only accommodate the 
unprecedented smallness of the light neutrino mass scale and all the peculiarities of the 
leptonic mixing pattern, but also provide specific new physics signals, thus admitting for 
a further experimental scrutiny. In some cases, such new physics effects could even be 
expected to be within the reach of near future experimental facilities. In this respect, the 
seesaw approach 5l-ll2l|. in which the smallness of the absolute neutrino mass scale is usually 
linked to a very specific type of high energy dynamics, represents a particularly plausible 
model-building paradigm. 

Each dynamical realization of the seesaw picture makes some kind of new physics effects 
appear, at least at a certain level. This, in turn, makes the neutrino sector an ideal probe to 
physics beyond the SM. For instance, the Majorana nature of the light neutrinos inherent to 
the seesaw framework provides characteristic lepton-number- violating signals at low energies 
like, e.g., the neutrinoless double beta decay, or, if kinematically accessible, same-sign di- 
lepton production at colliders, see, e.g., [13|. Similarly, besides neutrino oscillations, the 
flavour structure of the lepton sector can be tested in lepton-flavour- violating processes such 
as /i — 67 or, for example, trilepton collider events, c.f., fl^ . 

A full exploration of such new physics signals generally requires a very good knowledge 
of the leptonic flavour mixing angles governing the neutrino oscillation phenomena. In 
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particular, the smallest mixing angle 6*13, which plays a central role in the leptonic CP 
violation, is still to be determined (with a current 90% C.L. upper bound of sin^ 2^13 < 0.17 
reported by the CHOOZ collaboration 1^ and first indications of its non-zero value obtained 
by T2K jiel and MINOS (l7| experiments). 



The Daya Bay neutrino oscillation experiment [18[ is designed to perform a precision 
determination of ^13 with a potential to improve the CHOOZ limit by one order of magnitude. 
More specifically, the sin^ 26'13-sensitivity of Daya Bay is anticipated to reach 0.01 at 90% 
confidence level over the entire allowed range of Am^i, see Fig. 3.12. in 18|. At the assumed 
best fit point ^m'^-^ = 2.51 x 10~^eV^ the expected sensitivity is around 0.008 with 3 years 
of data. Such a highly ambitious goal relies on a very good control over the systematics, 
which is achieved by employing a unique set of eight identical detectors deployed at three 
different locations optimized for monitoring the antineutrino rates from the six reactors. A 
similar experimental setup is also adopted by the upcoming Double Chooz experiment |l^ 
and RENO 201] experiments. 

In combination with the large statistics due to the huge flux of antineutrinos produced in 
the nearby nuclear reactors, the unprecedented accuracy of this new generation of reactor 
experiments can make them sensitive to the new physics effects, at least at a certain level. 
For instance, if the new physics sector couples to hadrons and the relevant scale is not 
very high, one can expect non-standard interactions (NSI's) in the antineutrino production 
and detection processes as well as non-standard matter effects the antineutrinos experience 
throughout the propagation process. Similarly, new neutral fermions can mix with the three 
SM active neutrinos, which would result in an effective non-unitarity of the leptonic mixing 
matrix entering the relevant oscillations probabilities. 



The NSI's in reactor neutrino experiments have been discussed previously in, e.g., |2lL 12 2|. 
especially when the production and detection processes (and the corresponding non-standard 
effects) are assumed to be just inverse of each other. In particular, it has been shown that 
in such a case the NSI effects can be subsumed into a mere shift in the measured value of 
the effective mixing angle ^13. 

In this work, we study the NSI's in reactor antineutrino experiments in a general case 
when the assumption that the source and at the detector processes including the non- 
standard effects are just inverse of each other is dropped. This, in turn, leads to a specific 
distortion in the ratio of the antineutrino spectra measured in the far and in the near 
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detectors which can not be entirely transformed away by mere shifts in the relevant oscillation 
parameters, i.e., 6*13 and Amg^. We show that, under certain conditions, three years of Daya 
Bay running can be sufficient to provide a clear hint of non-standard neutrino physics. 

The remainder of the manuscript is organized as follows: In Sec. [TTl we present the 
general formalism and derive the relevant antineutrino survival probability formulas used in 
the subsequent analysis. Two basic scenarios corresponding to qualitatively different shapes 
of NSI's are specified in Sect. IIIIl and a detailed analysis of the observability of such effects 
at Daya Bay is performed in Sect. IIV[ Finally, we conclude in Sect. |Vl 



II. NON-STANDARD INTERACTIONS IN REACTOR ANTINEUTRINO OSCIL- 
LATIONS 

A. Non-standard interactions in the antineutrino sources and detectors 



In what follows, we adopt the standard SPD (source, propagation, and detector) ap- 



proach [23| to consider the antineutrino oscillation process in a reactor antineutrino exper- 
iment. In the presence of NSI's, the antineutrino states produced in the source as well as 
those observed in the detector can be treated as superpositions of pure orthonormal flavour 
states 

l3=e,fj,,T a=e,fM,T 

where the superscripts 's' and 'rf' denote the source and the detector, respectively. Note that 
there is no need to include the appropriate normalization factors in expressions ([1]) because 
we are going to be interested only in ratios of the survival probabilities in the near and far 
detectors where such factors cancel. 

The current experimental bounds on the NSI parameters mainly come from the lepton 
flavour violating decays ia — ?■ ^/37, the universality test of weak interactions and the invisible 
decay width of the Z-boson. Model-independent studies indicate that the upper limits on 



the NSI parameters e^'^ are typically in the ballpark of 10^^ to 10~^, see e.g. 2J] and 
references therein. 

To keep the analysis as general as possible, in what follows we shall consider two basic 
physical situations, namely. Eg = £^ (to be called case I, cf. Section IIII Al) and also Es 7^ £^ 
(case II, see Section IIII Bp and provide a detailed performance analysis of the Daya Bay 
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experimental setting in each case. 

Few comments are worth here: First, it is not very common to assume Eg 7^ as, in such 
a case, given the inverse microscopic nature of the relevant source and detection processes 
in standard reactor neutrino experiments, one usually concludes that these two quantities 
should be equal, see e.g. j2l|. This, however, is based on several implicit assumptions, 
in particular: 1) exact CPT (and Lorentz) invariance and 2) factorizability of the source, 
propagation and detection processes. As for the former, CPT violation (implying also the 
Lorentz invariance breakdown), for instance, makes it possible to have masses and mixings 
different for neutrinos and antineutrinos. Hence, the mixing matrix entering the detection 
process can be different from the one governing the antineutrino production and, thus, the 
standard oscillation formula is not valid even if all the NSI parameters are set to zero. 
Nevertheless, at least to the leading order in Eg and Ed , such an effect can be still accounted 
for within the standard SPD formalism by a suitable redefinition of the NSI parameters. 
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thus generating an effective departure from the assumed cf. formula (E 

Section FlI CI An interested reader can find further comments in, e.g., the recent work 25| 
and references therein. Needless to say, these considerations became especially relevant 
with the recent claims of observation of possibly superluminal neutrinos in the OPERA 
experiment 26|. Concerning 2), it is clear that as long as the antineutrino is produced 
in the source by the classical SM charged current interaction no departure from the basic 
case I setting can be expected because all the flavour-blind standard nuclear effects such 
as, e.g., the enhancement of scalar and/or tensor modes etc., simply factorize out and, 
hence, contribute only to the change of normalization of the total neutrino fluxes. In this 
respect, the Es 7^ case can be viewed as a model-independent parametrization of non- 
factorizable flavour-dependent new physics effects which, at low energies, can correspond to 
higher-dimensional terms in the effective lagrangian such as, e.g., {NNYj^V (which can be 
relevant also to other processes such as neutron star cooling, see, e.g., [27|) etc. Therefore, 
though possibly marginal from the perspective of the conventional extensions of the SM, we 
still find this setting worth a closer look, the more that Daya Bay can perform very well in 
this case, see Sections lIV D 21 and IIV D 3l 
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B. Non-standard interactions in the antineutrino propagation 



The propagation of antineutrino flavour eigenstates from the sources to the detectors is 
governed by the effective Hamiltonian 

H = Ho + H^ + H^si = ^f/MiagK, ml ml)U^ - diag(ycc, 0, 0) - Vccs"'* , (2) 

where e"^ is a Hermitian matrix parametrizing the NSI's throughout the antineutrino propa- 
gation and Vcc = V^GpNe arises due to effects of the coherent forward scattering in matter 
(with A'^e denoting the electron number density along the antineutrino trajectory). Bar- 
ring the irrelevant Majorana phases, the vacuum leptonic mixing matrix U is conveniently 



parametrized by three mixing angles Ou, 623 and 613 and one Dirac CP phase 6 28 1 



U 



"■^12^23 ~ C]^2'^23'^13^'^ C;^2C23 "~ '5l2'^23'^13^'^ ■^23'^13 
y '^12'^23 ~ '^12'^23'^13^"^ ~'^12'^23 ~ '^12'^23'^13^"^ '^23''13 J 



(3) 



with = cos and s^j = sinO^j (for ij = 12, 13 and 23). The full effective Hamiltonian 
(|2]) is diagonalized via a unitary transformation 

H = ^U*dmg {ml ml ml) , (4) 

where rhj {i = 1, 2, 3) denote the effective masses of neutrinos and U is the effective leptonic 
mixing matrix in matter. 

The size of the matter effect driven by the y/^GpNe term amounts to around 1.1 x 
10"^ eV^/MeV for Earth crust with density of 2.8 g/cm^. Even for the highest values of 
reactor antineutrino energies of around 10 MeV, this number is about 40 times smaller 
than the value of Aml/{2E) = 3.8 x 10"^ eV^/MeV and about 1100 times smaller than 
Aml2/{2E) = 1.2 X 10"^ eVVMeV. This indicates that Earth matter effects are very small 
and can be safely neglected. Hence, we take H Hq or, equivalently, set Vcc = in Eq. ([21). 

C. The antineutrino survival probability 

With the NSI effects at play, the electron antineutrino survival probability amplitude 
A{l4 ^1^-L) = Aee{L) is given by 

Aee{L) = (^|e-'^^|]^) = (1 + e^*)peA^p (1 + e^e, = [^ + + + ^"^^'lee ' 
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where L is the propagation distance and A is a coherent sum over the contributions of all 
the mass eigenstates z/j 

U^,U;,e-'^ . (6) 

i 

The antineutrino survival probability is then given by P{vl — ?■ yf) = \Aee{L)f. For com- 
pleteness, let us remark that a corresponding neutrino oscillation amplitude can be readily 
obtained from ([5]) with a substitution {U*,e*) — )■ {U,e). It should also be stressed that only 
the first row of 5^* and the first column of e'^ are relevant to the ee-type transition amplitude. 
Namely, the NSI parameters and e'^ involved in reactor neutrino experiment contain at 
least one flavour index e. 

Inserting formula ([6]) into Eq. ([5]) one arrives at the full antineutrino oscillation probability 

P{14 ^ ^) = ^ TJ^*-4: MJ'jn sin' (^J§^) +2 Yl MJ'Jn sin 

i,j i>j ^ ' i>j 

(7) 

where Am'^j = rrif — m|, and 

= UeiU^^ + ^ ] e^^U^iUi,^ + ^ ] E^gUeiU^i + ^ ] ^ e-y^ pe^^i^ pi ■ (^) 
7 7 7.P 

In the e'^'"' — )■ limit, Eq. ([7]) reduces to the standard survival probability. 

In this study, the quantity of our main interest is the third term in Eq. ([7]) which, being 
linear in the sine of L/E, does not play any role in the standard oscillation case. In this 
respect, a potential deviation from the "standard" quadratic-sine L/E dependence in an 
oscillation experiment can be interpreted as a hint of non-standard antineutrino interactions, 
in particular if such an anomaly exhibits the characteristic linear-sine L/E shape. 



Am^jL 
2E 



D. Series expansion of the antineutrino survival probability 

In practice, given the finite precision of the experimental inputs, it is very convenient 
to expand the survival probability ([7]) around the standard oscillation formula in terms of 
the relevant small parameters, in particular e^'^ which are all expected to be at most at the 



few per-cent level, c.f. 24 and references therein. In addition, ^13 is small compared to the 
other mixing angles (with the current CHOOZ upper limit of sin^ 2^13 < 0.17) and, hence, it 
amounts to another useful expansion parameter. Moreover, for the Daya Bay far detector. 



also the oscillation term Am|^L/(2i?) turns out to be at the level of 10~^ to 10~^ and, as 
such, it can also be viewed as a small quantity. 

Taking all this into account, we obtain the following expanded form of the relevant 
electron antineutrino survival probability 

P{Ul ^ ^) ~ P{ir, ^ I^)sM - 4 [Re {el^e-'' + ej.e'') + Re {el,e-'' + ei/') C23S13 

+ Re [el^el) si, + Re [elet) + Re {el/^, + e^^^,) .23C23] sin^ (^^) 
+ 2 [Im {el^e-'' + ej,e'^) .23^13 + Im (e^^e"'^ + e^/') c^ss,. 



+ 2 [Im {el^ + e'l^ C12S12C23 - Im [el^ + e^^) C12S12S23] 



2E 



(9) 



+ O 



^3 „3 _„2 / Aml-^^L 

^ 5'5i3,£ Si3,£:S^3,£:Si3 | 



2E 



2 / Am^^L 
''^M 2i? 



where P(z^e — ^ i^e)sM corresponds to the standard oscillation probability, i.e., the one without 
NSI's which is approximately given by 



1 - 4sf3 sm^ ' ''^ 



„3 2 
■=■13' ■'IS 



Am2iL 



2E 



10) 



Inspecting Eq. ([9]) one can recognize three qualitatively different non-standard contributions 
to -P(^'| — )■ )• ^he first two lines there is a CP-even term quadratic in sine of Am|2-Z^/(4£^) 
which, as expected, may affect the determination of the mixing angle ^13. The remaining 
three lines denote the CP-odd NSI effects corresponding to two different kinematical regimes 
characterized by Am|2 and Am2i, respectively. Notice that in the standard parametrization 
([3]), the Dirac CP-violating phase groups only with the former factor. It is also worth noticing 
that the term proportional to Am^i tends to be further suppressed in the "flavour-blind" 
setting (with = e^'^) because of the proximity of S23 and C23. 

It shall be noted that both the standard and the NSI transition probabilities depend on 
the neutrino mass hierarchy. For the standard oscillations, the hierarchy-sensitive terms are 
of the order of sl,Aml-^^L/{2E) and thus can be consistently neglected in Eq. ffTOj) . The 
NSI-dependent terms in Eq. (Q, however, contain a term linear in sine of Am\2L/{2E) 
which, indeed, differs in sign in the normal and in the inverted hierarchy schemes, respec- 
tively. Since, however, we do not expect any distinctive NSI features to be large enough to 
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discriminate among these two settings (although they would certainly differ in details), in 
what follows, we shall deliberately stick to the normal hierarchy case, i.e., assume Aml2 > 0. 



E. Notation and conventions 



In what follows we shall adopt the following parametrization: 




where the universal e index was dropped for simplicity. It is also convenient to define the 
source and detector phase averages $q, and differences A^^, respectively: 

<fa = i(0'+C), A0.^ 1(0^-0^). (12) 

The latter has a clear physical meaning: indeed, for all Acpa — ^ (together with — |£^|) 
one recovers a limit in which the non-standard antineutrino interactions in the detection 
process are of the same kind as those in the production. 



III. SPECIFIC SETTINGS 



In what follows, we shall discuss two simple but phenomenologically interesting shapes 
of NSI's and discuss the relevant effects in the reactor antineutrino experiments. 



A. Case I: e% = 

We start with the simplest case characterized by the assumption = e^* = le^le^''^" 
which corresponds to the situation where the production and the detection processes (in- 
cluding the associated non-standard interactions) are just inverse of each other. The relevant 
antineutrino survival probability (111) is then reduced to 

P(l4 ~ PiWe I^)sM - 4 {Saale;.!^ + cl^\er\^ + 2S2zC2'i\e^j\\er\ cos(0^ - (pr) 

+ 2si3 [szsk/,! cos(0^ -5) + caskri cos {(pr - 5)]} sin^ f \ . (13) 

Remarkably, the linear sine- dependent term in Eq. (jH]) vanishes and the NSI effects enter 
the survival probability as a mere global shift of the oscillation amplitude. This amounts to 
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FIG. 1: The oscillation probability for sin^ 26*13 = 0.1, 5 = and L = 1.8 km with no NSI's (solid 



line) and with NSI's corresponding to Case I in Sect. IIII Al with = e^* (dashed line) where we 
adopted = e,- = 0.02 (with both (f)^ and (f)^ fixed to zero). For the other oscillation parameters, 
the best-fit values have been assumed, c.f. Ref. ^. 



a shift in the "effective" reactor mixing angle 

■^13 = sl^+ sl^\e^\^ + cl^\er\^ + 2S2'iC2'i\e^\\er\cOs{(l)^- (t)r) 

+ 2si3 Nsk^tl cos(0^ -5) + caskri cos(0^ - 5)] . (14) 

Namely, the oscillation probability is given by the standard formula ( ITOj) with 6*13 replaced 
by the effective mixing angle ^13. Thus, there is no way to discriminate such an NSI effect 
from standard oscillations in reactor antineutrino experiments. It is also worth noting that 
the CP phase differences enter Eq. ( fT4l) via cosines only which is, indeed, justified by the 
CP properties of the survival probability in the setting under consideration. 

In Fig. [H we display the standard and the modified oscillation probability in the NSI 
presence as a function of the antineutrino energy in a detector at the "ideal" distance L = 1.8 
km (optimized for the highest count rate at ~ 4 MeV) from the source. The "depth" of the 
first oscillation minimum (the solid line for the standard oscillations) changes significantly if 
the NSI effects are turned on (dashed line); however, the energy of the minimum determined 
by the neutrino mass-squared differences remains essentially unchanged.^ 

^ Let us remark that the energies corresponding to the two relevant minima are not exactly the same due 
to the presence of the sub-leading terms proportional to (^Am^iL/AE)^ in Eqs. ([9|) and ([T0| . 
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Nevertheless, though the reactor antineutrino experiments in this case cannot distinguish 
the NSI's from a true mixing angle on their own, they can still provide a useful piece 
of information in combination with other types of experiments such as, e.g., accelerator 
experiments, superbeams, beta-beams, neutrino factories, etc. In particular, if these searches 
report different values of ^13, NSI's could be responsible for the mismatch. 

B. Case II: 4 ^ 

As already mentioned in Section Hi A[ we do not intend to confine ourselves entirely to the 
"canonical" case with = e^*, but rather prefer to keep the mind open also to the intriguing 

7^ ^a* possibility. Indeed, the distortion of the shape of the detected antineutrino spectra 
with respect to the standard oscillation picture expected with such a choice of the NSI 
parameters in the effective quantum-mechanical SPD picture can mimic (at the leading 
order) a wide class of new physics effects such net CPT violation or, for instance, 

non-factorizable beyond-Standard-Model effects in the production processes, cf. Section [ITl 

As a consequence, the terms linear in sine in formula ([9]) are exposed and the relevant 
NSI effects can no longer be completely subsumed into a shift of the effective mixing angle 
^13. This, besides the change of the "depth" of the first oscillation minimum (cf. Figured]), 
leads also to a shift in its energy, as illustrated in Figure [2J In particular, the dip can be 
shifted by as much as one MeV in both directions, depending on the specific choice of the 
NSI parameters. 

In what follows we shall focus on two specific realizations of this setting, namely, the case 
when the magnitude of the NSI parameters differs between the production and detection 
processes (Case Ila) and the case when the relevant NSI parameters are of the same size but 
differ by their phases (Case lib). Both these cases are studied numerically in Sect. IIVI 

1. Case Ila: Non-standard interactions in source only 

Let us exemplify the first option on a specific setting where the NSI's exhibit themselves 
only in the production processes, i.e., taking e'^ = 0. Given that, the general formula (jH]) 
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0.85 L 



FIG. 2: The theoretical oscillation probability with no NSI (solid line) and with NSI's (dashed 
and dotted lines) for L = 1.8 km. For the sake of illustration, we adopt a flavour-universal scheme 
with all the relevant NSI parameters at the same level of \£^,t\ = 0.04 with = ^. The other 

neutrino mixing parameters are the same as in Fig. [TJ The dashed line corresponds to $^,r = f 
whereas the dotted one to <^^^t- = |7r, respectively. 



simplifies into 

- 4si3 [sasl^^l cos(</)^ - 5) + CsskrI cos(0^ - 5)] sin^ 

+ 2si3 [■523^^1 sin(0^ - (5) + C22.\er\ sin(0^ - 5)] sin 

+ 2si2Ci2 (c23|e^| sin0^ - S23kr| sin0^) f — j , (15) 

where we used Sa = and 0^ = 0a- It is worth noting that the genuine NSI effect (due 
to the last two terms) is proportional to sines of differences of the Dirac CP phase 5 and 
the CP phases of the NSI parameters 0^, as expected for a CP- violating effect beyond the 
standard oscillation picture. 



f AmLL\ 



f AmLL\ 



2. Case lib: Same-size source and detector effects with different phases: |e^| = |e^|, A0q 

An interesting complementary setting is obtained if the magnitude of the source and 
detector effects are assumed to be equal so that the NSI effects can only be distinguished due 
to the mismatch between the corresponding CP phases 0^ and 0^. With = \e^\ = \ea\, 
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the general formula ([9]) receives a rather symmetric form^ 

— 4 {S23|e^|^ COS 2A{/)^ + C23|£:rP COS 2A(/)t- 
+ 2C23S23k/.|kr| COs(A0^ + A</)^) COs(<l>^ - 

f Am? L\ 

+ 2si3 [s2z\e^,\ cosA0^cos(<l>^ -5) + C2z\er\ cosA0^cos($^ - 5)]}sin^ ( \ 

— 2 {s23|£:^P sin2A0^ + C23|er|^ sin2A0T 

+ 2023523^^.1 kr I sin(A0^ + A0^) cos($/, - <l>^) 

+ 2si3 [s23k^| sin A</)^ cos($^ - 5) + C23|£:^| sin A</)^cos($^ - 5)]} sin 

— 4si2Ci2 (c23k^| sin A0^,cos<l>^, - S23|£:r| sin A0^COS$^) 



Am\2L 
2E 

AmhL\ , , 



where the notation specified in Eq. ( |T2l) has been used. Again, the relevant phase differences 
in the genuine NSI terms enter in sines. Furthermore, the formula above can be simplified 
to a yet more compact form 



^ + 2A0,)+4M^sm( J 

1e 



4^ 4k/.rsin [^-^ + 2A0^j + ^k.psin ( + 2A0. 

+ 2C23S23k/.|kr| slu g + A0^ + A(/)^ COs(<l>^ - 



(Am} L \ 
+ 2gi3g23|gM|sin ( + A0^j cos($^-^) 

, , ( Am%L ^ \ , ^,-\ ( Am\^L\ 
+ 2si3C23kr| sin ( + A0, 1 cos(<l>, - 5) I sin ( 1 

f Am^ L\ 

- 4si2Ci2 (023^^1 sin A0^ cos $^ - S23kT I sin A0^ cos f j ,(17) 

which does expose the "kinematic" role of the phase differences A0q, and the "amplitude 
modulation" role of their averages 

Let us also remark that for Acf)^ — > (when the symmetric setting with e'^ = e'^'^ is 
recovered) the last term vanishes and, as expected, the other terms conspire to yield a mere 
shift in the effective mixing angle ^13 identical to that given in formula (fT^ . This provides 



^ Note that the coefficient of the last term in Eq. ((T5)) is optically different from the same in Eq. ^ which 
is due to the utilized goniometric identity for a difference of two sines and the definition of A(/)q,. 
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a nice consistency check of the results. A simple numerical analysis of both Case-IIa and 
Case-IIb settings is given in Sect. [IVl 



IV. NON-STANDARD ANTINEUTRINO INTERACTIONS AT DAYA BAY 
A. Experimental setting 

The Daya Bay neutrino experiment 



18| is designed to perform a precision measurement 
of using antineutrinos produced by the reactors of the Daya Bay Nuclear Power Plant 
(NPP) and the Ling Ao NPP. In the detectors, antineutrinos from the reactors are captured 
via the inverse beta-decay process, and the deficit from the expected 1 /L^ dependence is to 
be interpreted as a signature of neutrino oscillations. In particular, near and far detectors 
are employed in order to suppress the systematic uncertainties related to the antineutrino 
flux from the reactors. 

The Daya Bay measurement of sin^ 26'i3 is expected to reach the sensitivity of the order of 
0.01, an order of magnitude better than the current CHOOZ limit sin^2^i3 < 0.17. Besides 
a high-quality determination of the relevant standard neutrino oscillation parameters, Daya 
Bay can be rather efficient in improving some of the current constraints on physics beyond 
the SM. 

In order to estimate the NSI effects possibly observable at Daya Bay, we perform a basic 
numerical analysis making use of a simple model of the detected neutrino spectra. There are 
three pairs of nuclear reactor cores of a total thermal power of 17.4 GW at the experiment 
site, namely, Daya Bay (DYB), Ling Ao (LA) and Ling Ao II (LAII), providing electron 
antineutrinos to three detectors, two near ones called Daya Bay (DYB) and Ling Ao (LA) 
with 40 tons and a far detector (FAR) with 80 tons of a Gadolinium-doped liquid scintillator, 
respectively. A more detailed breakdown of the relevant Daya Bay parameters can be found 
in TABLE E 

B. A simple model of the spectra 

For the sake of simplicity, we shall assume that each pair of neighboring cores constitute 
a single point source. The average energy release per one fission Ep is anticipated to be 
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DYB layout 


geometry 


expected daily event rates 


Detectors\Cores 


DYB 2x2.9 GW 


LA 2x2.9 GW 


LAII 2x2.9 GW 


mods.xDR [18] 


simulated 


DYB (40 t) 


363 


1347 


1985 


2 X 930 


2 X 890 


LA (40 t) 


857 


481 


1618 


2 X 760 


2 X 790 


FAR (80 t) 


1307 


526 


1613 


4 X 90 


4 X 90 



TABLE I: The basic Daya Bay experimental layout [1^ and expected daily antineutrino event 
rates: distances in meters between detectors (in rows) and centers of pairs of the neighboring 
reactor cores (in data columns 1-3) and numbers of anticipated antineutrino events per day at each 
site (in data columns 4 and 5). In column 4 we display the data quoted in the DYB proposal (cf. 
Table 3.8. in jlj]) multiplied by the number of modules, the numbers in column 5 correspond to 
our numerical analysis described in Section IIVI Indeed, the simulated event rates in all cases lay 
within 4.5% of the nominal Daya Bay values, thus justifying the relevance of the simplified model 
of the antineutrino spectrum as well as the expected statistical error levels. 



around 200 MeV so the estimated number of fissions per second in each reactor site Np 
is 

Np = 2 Pr/Ef = 2 X 2.9 GW/200 MeV = 1.8 x 10^° s"^ , (18) 

where the extra factor 2 counts the number of reactor cores per site and Pt stands for the 
thermal power of each core. For the spectrum of the antineutrino flux per fission we shall 
use the approximate formula given in Ref. 29|] (for E in MeV): 



^ = exp (0.87 - 0.16E - 0.091^^) MeV"^ 



(19) 



Antineutrinos interact with the free protons in the scintillator via the inverse beta decay 
process + p n + . The cross-section of this reaction has been calculated in Ref. jSOj 
to be 



(j{E) = 9.52 X 10 



-48 



{E - {rrin - rrip)) \J {E - (m„ - mp)f - ml MeV" 



m 



(20) 



with the energy threshold Eq = 1.8 MeV. There are 6.29 x 10^^ free protons in a cm^ of the 
scintillator of density p = 0.86 g/cm^ 3l|. Therefore the number of targets per one ton of 
the scintillator is Nt = 7.3 x 10^*^ ton~^ 
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dN 

— [10^ MeV"'] 
dE 

6 




E(MeV) 



FIG. 3: The expected shape of the detected antineutrino spectrum in the DYB (near) detector 
after three years of running without NSI's. 

The antineutrino survival probability Piiy^. — )■ u^) is a function of energy, propagation 
distance, oscillation parameters and, in general, also the NSI parameters. The expected 
total number of antineutrino events in the detector D (with D =DYB,LA,FAR) with mass 
Md after three years of running can be estimated as 



dN, 



D 



dE 



tx NtX MdX Ceff X Npx <^{E) xa{E)x ^ ^ ^^"^ 



_R=DYB,LA,LAII 



(21) 



where t = 3 x 365 x 24 x 3600 s is the duration of a three-years' run. We sum over three 
reactor sites and use Lor for the distance between the detector D and the reactor site R, 
c.f. TABLE [H In addition, we adopt a detection efficiency coefficient Ceff = 0.78 {isl. As 
an example, we depict in Figure [3] the expected spectrum of antineutrinos detected in the 
DYB detector. It is worth noting that the highest event rate corresponds to ^ 4 MeV. 

The quantity of our main interest is the ratio of the antineutrino energy spectra between 
the considered far and near detectors, which can be obtained readily from Eq. (I2ip . It is 
expected that, due to a similar design of the far and near detectors, the systematic uncer- 
tainties associated to, e.g., the absolute flux determination, can be greatly reduced in the 
ratio of the energy spectra. However, in order to fully account for all the systematic uncer- 
tainties, e.g., the backgrounds, energy miscalibration, detection efficiencies etc., a complex 
simulation of the Daya Bay experiment is necessary. This, however, is a formidable task in 
general, so in what follows we shall consider mainly the statistical uncertainties (Sections 
IIVCI and IIVD|) and only later on (in Section lIVEp we demonstrate that the changes due to 
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FIG. 4: Ratio between the FAR and the DYB detectors for standard oscillations with statistical 
errors only. The dotted line corresponds to the case of sin^ 2^13 = 0. The neutrino mixing 
parameters are the same as those used in Figure [TJ 



the (leading-order) systematical effects do not inflict any signiflcant changes on these results. 



C. The analysis 

To assess the observability of NSI's at Daya Bay in practice, we perform a simple numer- 
ical analysis along the following lines: we choose 15 energy bins from 1.8 MeV to 8 MeV 
in order to have approximately the same statistics in all bins which are 1-4 times wider than 



the energy resolution 15%/ \/ -E(MeV) 18|. In each bin, we use Eq. f l2T]) to calculate the 
ratio R of the antineutrino energy spectra between the far and near detectors (for sake of 
illustration, from now on we shall focus in particular onto the FAR and the DYB detectors). 
In the case of the standard neutrino oscillations, the expected shape of this ratio between 
the FAR and the DYB detectors is depicted in Figure HI 

Since the uncertainties of 6'23, 6*12 and Am^i are not expected to play any signiflcant 
role in the ratio of our interest, we shall flx these parameters to their central values. This 
is not the case of Am32 because the uncertainty in this parameter (quantifled by cr^mlj) 
mimics the effects of the NSI's, namely, it also shifts the position of the flrst minimum in 
R. However, with the increasing precision of the Am^g determination, these effects become 
less important. Therefore, in what follows, we shall mainly focus on two speciflc situation 
corresponding to different choices of (yi^m\^- iii case we take ^Aml^ = 0-09 x 10~'^eV^ ^ 
as the current experimental value while in the other "ideal case", we push o"Am§2 down to 
0.025 X 10~^eV^, respectively (which can be viewed as an optimistic expectation for the 
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uncertainty in the atmospheric mass-squared difference in several years from now) . 

For each specific choice of the relevant NSI parameters there are only two unknown 
parameters left in R, namely, S13 = sin6'i3 and Amgg, c.f., Eq. (ITU]) . Denoting the i-th 
bin value of R (as a function of S13, Amgg, and e*^) by -Rj(si3, Amgg, e**, e*^), we attempt 
to fit the simulated data by the NSI null-hypothesis corresponding to the case when Ri 
is calculated for standard oscillations with some effective values of the relevant oscillation 
parameters, namely, R^ (§13, Arrigg). This is done by minimization of the function 



15 
i=l 



R, (sl3,Ami2,£^£'^) -/?°(Si3,A 



™2 



m 



32) 



/ Am32 — Am32 



0"data 



(22) 



with respect to §13 and Am32, where (Xdata denotes the three- years' run statistical error (s). 
The §13 parameter has been left free (to be determined by Daya Bay) whereas Am^g should 
obey the existing experimental constraints. The value of iii the minimum (Xmin) then 
quantifies the likelihood that the Daya Bay data could be fitted by the standard oscillation 
formula. For two fitted parameters, xLin — 2-3 and 4.61 corresponds to 68% and 90% C.L., 
respectively. 



D. Results 



For the sake of simplicity, in what follows we shall consider only the "flavourless" versions 
of the oscillation probability formulas relevant to the three cases of our interest, namely 
Eqs. (fT3l) . (fTS!) and (fTTT) . This amounts to setting 0^ = 0^ = and |5^| = |£:^| = \e\ 
everywhere. Let us recall that, besides the standard oscillation parameters, in Case I and 
Case Ila the relevant input NSI parameters are, namely, the (universal) magnitude of the NSI 
effects l^l and the corresponding CP phase (more precisely, the phase difference 0' = — 5 
where S denotes the leptonic Dirac CP phase) while in Case lib the NSI parameters entering 
the survival probability are $ (again, it is rather $' = $ — 5) and A0. 

Let us also reiterate that only the statistical errors have been taken into account in the 
current analysis. A complete study including also the systematic uncertainties would require 
a complex simulation of the Daya Bay experiment. This, however, is beyond the scope of 
the present study. 



18 



R 



R 



E(MeV) 




E(MeV) 



0.215 



-2^ 3 4 5 6 7 



8 



0.215 



8 



0.205 



0.205 



0.195 



0.195 



0.185 



0.185 



FIG. 5: Sample fits of the ratio between the detected antineutrino spectra in the FAR and DYB 
near detectors in the case that the production and detection process (including the relevant non- 
standard interactions) are just inverse of each other (Case I, Sect. IIII Ap . i.e., e'^ = e"^^ . In the left 



\e\ = 0.02 and = instead. The dotted lines correspond to the standard oscillations without 
NSI's, while the solid lines are the fits based on the standard oscillation survival probability (jlOp 



sin2 2^i3 = 0.138 (right panel). 
1. Case I 

As we argued in Sect. IIII A[ in the symmetric setting with = e'^'^ , the NSI effects cannot 
be distinguished from the pure standard oscillations. Even if the underlying mixing angle 6*13 
is zero, one can still fit the data with a standard oscillation curve corresponding to a nonzero 
value of the effective mixing angle ^13 given by formula (fT4l) . A pair of representative plots 
depicting the expected ratio between the FAR and the DYB antineutrino spectra in this 
situation are given in Figure O One can see that the data are well fitted by the standard 
oscillation formula with just the effective mixing angles different from their "true" values. 



panel, sin^ 2^13 = 0, |e| = 0.04 and (f>' = have been used; in the right panel, sin^ 2^13 = 0.1 



used in Eq. (j2ip with the effective mixing angles given by sin^ 2^13 



0.013 (left panel) and 
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2. Case II a 



In the "flavourless" setting, the relevant Case-IIa formula (fT5|) simplifles into 




( 



AE 



) 



(23) 



where cj)' = (p — 5. Note that in most cases the last term can be neglected due to the 
experimental proximity of ^23 to ^. 

We also stress that for 0' — )■ or vr the leading NSI contribution corresponding to the 
sine-squared term above essentially mimics the effects of standard oscillations with a shifted 
mixing angle because, besides the last negligible term, there is no net NSI induced CP- 
violating effect. One can see this on the left panel in Figure [6] where, indeed, the data can 
be fltted by the standard oscillation formula with just a shifted effective mixing angle ^13. 
However, the change is still proportional to S13 and thus no shift is induced if the underlying 
6'i3 happens to be zero^. 

However, for non-trivial 0', the NSI effects can no longer be subsumed into a pure shift 
in ^13 and the standard oscillation formula no longer flts the data even if one admits for 
a certain variation in Amgg, see the right panel in Figure [6l Thus, in this case, one can 
in principle attempt to constrain the \e\ and 0' parameters, at least in some parts of their 
parameter space. 

In Figure [7] we present the relevant exclusion regions for these parameters. Therein, one 
can observe an interesting vr-periodicity in 0', which can be understood from the shape of 
the second correction in formula fl23|l . Remarkably enough, even with variable Am32 and 
l^l as low as 0.02, in some cases the NSI effects can be distinguished from the standard 
oscillation at 90% C.L. 

^ As we shall see, this is different from the Case lib setting studied in Sect. IIVD 31 where a nonzero value 
of the effective mixing angle can be generated even for ^13 = 0. 
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FIG. 6: Sample fits of the ratio between the detected antineutrino spectra in the FAR and DYB 
near detectors in the asymmetric setting where the NSFs are assumed to affect only the production 
process (Case Ila, Sect. [IIIB ip . We adopt sin^ 2^13 = 0.1 and |e| = 0.04 in both panels. Further- 
more, in the left panel, 0' = is assumed, while in the right panel we take cp' = ^ (maximal CP 
phase difference). The dotted lines correspond to the standard oscillations without NSI's, while 
the dashed lines show the fitted curves with the effective mixing angle sin^ 2^13 = 0.135 (left panel) 
and sin^ 2^13 = 0.105 (right panel). In addition, the solid lines stand for the fitted curves with two 
parameters ^13 and Arh^2- the left panel, the solid line coincides with the dashed line, whereas 
in the right panel it does not and the best fit corresponding to the values sin^ 2^13 = 0.109 and 
Am|2 = 2.20 X 10"3 eV^ requires a significant shift in with respect to its central value. 

3. Case lib 

In the more general case when both the source and detector effects are present there is 
an extra set of parameters at play associated to e^, i.e., the relevant magnitudes \e'^\ and 
also the extra detector NSI phases which combine with the source ones into the phase 
averages <I>q, and the phase differences A0q,, c.f. Eq. ( |T2l) . As before, we will assume a 
"flavourless" form of NSI's and deliberately put \e'^\ = \e^\ = \e\ in order to simplify the 
numerical analysis. Then formula ( 1T6|) reduces to 



- 4(s23 + C23)^|£:|^sin 



sm 



+ 2A(p 



- 8Si3(s23 + C23)k|cOS$'sin ' * ' 



AE 



sm 



AE 



4Si2Ci2(C23 - S23)|£:| sinA0cos<I> ( 



(24) 
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FIG. 7: Regions of parameters (on the right hand side of the curves), where the Daya Bay exper- 
iment can disfavour the standard oscillation hypothesis at 68% (thin curves) and at 90% (thick 
curves) C.L. for sin^ 2^13 = 0.1. The solid curves stand for a standard two-parameter fit with 
variable ^13 and ^^mj^ — 0-09 ^ lO^^eV^ while the dashed curves correspond to the fit 

with o"Am§2 pushed down to 0.025 x 10~^eV^. 

where, again, $' = $ — 5. As before, the last term is negligible for 623 close to 7r/4. 

For the sake of illustration, in Figure E] we show two specific examples of the i?-fits 
obtained in Case lib. There, the data are fitted by the standard oscillations, first with 
variable 613 and (solid lines) and then also with only 6'i3 as a free parameter (dashed 

lines) . 

In Figure [9l the exclusion plots for the A0 and \e\ parameters are given for sin^ 26*13 and 
two specific choices of The sensitivity in \e\ is similar to that observed in Figure IH] for 
Case Ila. Notice, however, that the two leading corrections in Eq. ( [T6l) have a very different 
A0-periodicity. The former is vr-periodic in A0 while the latter is ^-periodic in A(f). The 
reason is easily seen from the analytic shape of the relevant survival probability (!24|) . Indeed, 
for $' = 0, the second correction in formula flMl) dominates over the first one while it is the 
other way round for $' = ^. 

Regions in the $' — A0 plane where Daya Bay experiment could distinguish non-standard 
effects from standard oscillations (at 68% and 90% C.L.) are shown in Figure[IO]for different 
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FIG. 8: Sample fits of the ratio between the detected antineutrino spectra in the FAR and DYB 
near detectors in Case lib. We adopt sin^ 2^13 = 0.1 and |e| = 0.02 as well as Ac/) = | in both 

is assumed, while in the right panel we put = |7r. 



TT ; 
2 



panels. Furthermore, in the left panel, 
The dotted lines correspond to the standard oscillations without NSI's, while the dashed lines show 
the fitted curves with the effective mixing angle sin^ 2^13 = 0.103 (left panel) and sin^ 2^13 = 0.092 



(right panel). In addition, the solid lines stand for the fitted curves with two parameters, i.e., 
sin2 2^i3 = 0.105 and Am|2 = 2.20 x 10"^ eV^ in the left panel, and sin^ 26*13 = 0.094 and 
Am§2 = 2.72 X 10"^ eV^ in the right panel. 



values of \e\ and sin^ 26*13. If the value of the underlying 6*13 is close to the CHOOZ limit 
(sin^ 2^13 < 0.17) and |£:| = 0.02 then the region is relatively large, see the upper-left panel in 
Figure [TOl With decreasing |5| (from left to right) or 6*13 (from up to down), the observability 
domains become naturally smaller. 

The possible NSI effects in an independent Daya Bay determination of the standard 
oscillation parameters 6*13 and Am^^ illustrated in Figure [11] One can see that, at least 
in some cases, the corresponding global best fit point can differ significantly from the "true" 
values of these parameters, potentially leading to a tension between Daya Bay and other 
experiments. 

Yet another comment is in order here. As we have seen in Case Ila (c.f.. Sect. IlVD 21) . 
with source effects only there is no way to end up with a significant effective ^13 if the 
underlying ^13 was zero, while here one still gets ^13 7^ even for ^13 = due to the first 
term in Eq. flM|) . Such a qualitative difference in the behavior of these two settings can be 
heuristically understood as follows: In the former case, there are effectively only two small 
parameters (with their corresponding CP phases) at play, namely \e'^\ and S13 while there 
are three such quantities in the latter case, in particular \e^\, \e'^\ and S13. In Case Ila, 
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FIG. 9: Regions of parameters (on the right hand side of the curves), where the Daya Bay exper- 
iment can disfavour the standard oscillation hypothesis at 68% (thin curves) and at 90% (thick 
curves) C.L. for sin^ 2^13 = 0.1. Here we assume = (left panel) and = ^ (right panel). 
As before, the solid curves stand for a standard two-parameter fit with variable ^13 and Am|2 for 
•^Am^j ~ ^'^^ ^ 10~^eV^ while the dashed curves correspond to the fit with CAm^j Pushed down to 
0.025 X lO-^eV^. 

there is thus only a single relevant phase difference governing the CP-even effects [due to 
the first correction in Eq. (125]) ] which, however, becomes ill defined in the S13 — )■ limit, 
and thus its effect can be "rotated away" . Remarkably, this is not so in Case lib since there 
is an observable phase difference A0 left even in the S13 — )■ limit and the corresponding 
contribution to the effective 613 due to the first term in Eq. ( 124|) cannot be transformed out. 



E. Effects of leading-order systematics 

Finally, let us argue that the leading systematical effects do not change the results ob- 
tained in the previous sections in any significant way. 

The main sources of systematical uncertainties in the Daya Bay setting are related to the 
reactors (power, spent fuel, location), detectors (energy miscallibration, target mass, detec- 



tor efficiency) and, of course, backgrounds (accidental signals, ^He ^Li, fast neutrons) |18 |. 
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Taking the full advantage of the "near+far" detector setting one can approximate the 
leading-order systematic uncertainties (namely, the neutrino flux uncertainty and the un- 
certainty in the detector masses) as a relative change of the measured far-to-near ratio of 
the detected antineutrino energy spectra R, cf. Section llVCi For the sake of simplicity 
we shall assume that this change (to be denoted by K), as well as its uncertainty ax, are 



we 



energy-independent at the leading order. Following the detailed discussion given in 
shall adopt a conservative value of ax = 0.6% for the calculation. 

The argument above makes it possible to implement the leading-order systematics by 
simply extending the original formula fl2^ into 
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^data 



(25) 



+ 



Am32 — Am|2 \ f K — 1^ 



(^Aml^ I \ 0-K 



and marginalizing over S13, AmL and K. With three degrees of freedom at play the 68% 



CL. and 90% CL. values now correspond to = 3.53 and = 6.25, respectively. 

Let us illustrate the smallness of the changes inflicted by the variation of K on, e.g., 
the situation of Case lib studied in detail in Section IIVD 31 In Figure [12] we demonstrate 
the shift in the solid contours displayed previously in Figure M (where these were obtained 
for (^Am^j = 0.09 X 10~^eV^ with only statistical uncertainties taken into account) due to 
the systematic effects. The dashed lines in Figure [12] demonstrate the slightly reduced 
sensitivity of the Daya Bay if the leading systematics is taken into account. To conclude, 
the systematics does not hinder the Daya Bay's sensitivity of to the new physics effects and 
the discovery reach remains safely statistics-dominated. 



V. SUMMARY AND OUTLOOK 



In this work, we have performed a detailed analysis of the non-standard antineutrino 
interaction effects in the Daya Bay short-baseline reactor antineutrino experiment. 

The NSI's in reactor antineutrino experiments can exhibit themselves in various ways 
depending on the character of the underlying physics. If, for instance, the non-standard 
interactions in the production and detection processes happen to be exactly the same, i.e., 

_ ^df^ j^g^ effect consists in a shift in the depth of the oscillation dip in the measured 
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ratio of the far and near detector antineutrino spectra corresponding to the extracted value 
of the mixing angle 613. Thus, in this case, the NSI effects can not be distinguished from 
the standard oscillations 22 1. 



If, however, this assumption is relaxed, owing to, e.g., non-standard multi-body interac- 
tions in the source, the measured antineutrino spectra are distorted in a specific way and 
become incompatible with the standard oscillation interpretation - besides the change of 
the depth of the first dip, also its energy position is shifted. This can be only partially 
accounted for by the standard oscillation formula if the extracted values of the mixing angle 
^13 and, in particular, the corresponding mass-squared difference Amig, are both allowed 
to differ significantly from their genuine values. However, in practice, the effect can not be 
entirely subsumed into a shift in the 6*13 — Amgg plane due to the strict constraints on these 
parameters from other measurements. 

In Sect, ini we have derived general formulas for the oscillation probabilities including the 
non-standard effects in the antineutrino production and detection processes, arguing that 
the matter effects throughout the antineutrino propagation do not play any significant role 
in short baseline reactor neutrino experiments such as Daya Bay. 

In Sect, mil we specified the setting of our main interest corresponding to three different 
configurations of the NSI parameters. In Sect. IIVI we performed an illustrative numerical 
analysis of these settings based on an empirical model of the reactor antineutrino spectrum at 
Daya Bay assuming for simplicity that the NSI effects are flavour blind. Taking into account 
the statistical uncertainties corresponding to three years of running, we have studied how 
the NSI's could modify the antineutrino energy spectra and the measured values of the 
neutrino mixing parameters in practice. We observe that, under certain conditions, the 
Daya Bay experiment can provide hints of such non-standard effects at more than 90% C.L. 
The leading-order systematics has been discussed in brief and it has been shown that it does 
not play any decisive role in the expected Daya Bay new-physics sensitivity. 

We should also stress the important complementary role the long baseline experiments, 
such as, e.g., accelerator experiments, superbeams, beta-beams or a neutrino factory, could 
play. Namely, if 61^ or Am|2 as determined by Daya Bay differ significantly from the other 
results, one would have to take the NSI effects seriously as one of the possible sources of such 
a discrepancy. In that case, a combined analysis of the Daya Bay data together with the 
data from the other experiments, including the NSI's of the kind considered in this study. 
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would be of utmost importance. 
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FIG. 10: Regions of parameters (interior), where the Daya Bay experiment can disfavour the 
standard oscillation hypothesis at 68% (thin curves) and at 90% (thick curves) C.L. for \e\ = 0.02 
(left column) and \e\ = 0.01 (right column) and sin^ 20i3 = 0.1 (upper row) and sin^ 20i3 = 0.05 
(lower row). Again, the solid curves stand for a standard two-parameter fit with variable ^13 and 
Am^2 ^^ml^ — 0-09 ^ 10~^eV^ while the dashed curves correspond to the fit with fAm^j Pushed 
down to 0.025 x lO-^eV^. 
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FIG. 11: The effects of the non-standard interactions in the determination of the standard os- 
cihation parameters ^13 and Daya Bay after 3 years of running. The upper cross 

denotes the assumed "true" values of the standard oscillation parameters sin^ 2^13 = 0.1, 
— 2-45 X 10~'^eV^. Turning on the NSI parameters (fixing, for instance, |e| = 0.02, 
$ = 'it/2, A(f) = 7r/2 in case lib, c.f., Sect. IIIIB 2^ . the best standard oscillation fit is shifted to 
sin^ 2^13 = 0.105 and Am32 = 2.20 x 10~^eV^ (the lower cross) and the corresponding Xmin ~ ^^-^ 
indicates a significant incompatibility between the Daya Bay data and the standard oscillation hy- 
pothesis. We display three solid curves depicting the levels around the best-fit point; from thick 
to thin, = 20, 40 and 60, respectively. The shaded bands depict the pull (due to the second 
term in Eq. (j22p l inflicted by A.rh^2 departing from the "true value"; the dark/light and light /white 
boundaries enclose the Am^j = (2.45±0.09) x lO^^eV^ and Am^g = (2.45±0.18) x lO^^eV^ regions 
(about la and 2a), respectively. 
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FIG. 12: The Daya Bay new-physics sensitivity with both statistics and systematics taken into 
account (in dashed Unes) as compared to the results obtained previously with solely statistical 
uncertainties (in solid lines), cf. Figure [9l As usual, the thick and thin lines correspond to the 
90% C.L. and 68% C.L. contours, respectively. All the NSI and oscillation parameters are fixed to 
the values employed in Figure [9l i.e., sin^ 2^13 = 0.1, = (left panel) and = ^ (right panel), 
with a^^2 = 0.09 x IQ-^ eV^. 
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